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THE MORPHOLOGY OF A MULTI-BUBBLE SYSTEM 
IN THE IONOSPHERE 


I. INTRODUCTION 

The behavior of Ionospheric plasma has been extensively studied with 
the aim of understanding various ionospheric phenomena including equatorial 
spread F (ESF) and plasma cloud strlations (see the following reviews: 
Ossakow, 1979; Fejer and Kelley, 1980; Ossakow , 1981; Kelley and McClure , 
1981; Ossakow et al. , 1982). In particular, ESF is thought to be initiated 
by the Rayleigh-Taylor instability, first proposed by Dungey (1956). In 
the context of this idea, the plasma density depletions resulting from the 
instability acquire upward polarization induced E_ x B_drift velocities. A 
large body of literature has since been developed to describe the linear 
and nonlinear properties of Rayleigh-Taylor plasma density irregularities 
under various assumptions [ Haerendel, 1975; paisley et al. , 1972; 
Chaturvedi and Kaw , 1975a,b; Hudson and Kennel, 1975.]. In particular, 
considerable attention has been given to the morphology and motion 
[ Scannapleco and Ossakow, 1976; Ossakow and Chaturvedi, 1978; Ott, 1978; 
Hudson , 1978; Anderson and Haerendel , 1979; Ossakow et al. , 1979; Zalesak 
and Ossakow, 1980; Zalesak et al., 1982] of plasma depletions ("bubbles") 
in the equatorial ionosphere. Moreover, a number of observations [ Woodman 
and LaHoz, 1976; Kelley et al., 1976; Hanson and Sanatani, 1971; McClure et 
al. , 1977; Szuszczewicz et al. , 1980, 1981] have indicated the presence of 
rising plasma bubbles. 

Another phenomenon of Interest is that of the strlations in 
(artificial) plasma "clouds" (density enhancements). This effect has been 
attributed to the x gradient drift instability [ Linson and Workman , 
1970; Volk and Haerendel, 1971] and appears to be amenable to treatments 
similar to ESF [Scannapleco and Ossakow , 1976; Scannapleco et al., 1976; 
Ossakow and Chaturvedi, 1978]. The Rayleigh-Taylor instability and the E_x 
B_ gradient drift Instability are both interchange modes that may occur in 
the leading edge of plasma bubbles and backside of clouds respectively. In 
both cases, the resulting density depletions and enhancements are thought 
to drift by the polarization induced E_ x B_ drift. The electric field is 
produced by the polarization of the plasma across the earth's magnetic 
field. The essential ingredient Is the small but nonzero ion-neutral 
collision frequency. 

Manuscript approved November 24,1982. 
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As a result of the Initial Instability such as the Rayleigh-Taylor 
Instability, an array of density depletions are formed with the wave "“ctor 
perpendicular to the earth's magnetic field. Moreover, the leading edge of 
a bubble itself is Raylelgh-Taylor unstable resulting in further 
"bifurcation" [ Ossakow and Chaturvedl, 1978]. Similarly, the x _B 
Instability causes the backside of an initial cloud to striate, forming 
"flnger-llke" structures. Thus, the plasma depletions and enhancements 
typically occur In multitudes. Indeed, McClure et al., (1977) and 
Szuazczewlcz et al« , (1980;1981) have given evidence for multiple bubbles. 
In the previous theoretical and numerical simulation works [ Ossakow and 
Chaturvedl, 1978; Anderson and Haerendel , 1979; Zalesak and Ossakow , 1980; 
Overman et al., 1982], the morphology and evolution have been studied using 
one-bubble models or a uniformly distributed array of bubbles (clouds). As 
a result, the intrinsic influence of neighboring bubbles on each other has 
not been well quantified. Because the bubbles rise due to the polarization 
Induced E_ x B_ drift, the electric field configuration In and around the 
bubbles is of central importance. 

In the present paper, we study the structure of electrically 
Interacting multi-bubble systems and seek to Identify the nature and 
effects of the mutual interaction. For this purpose, we start with a 
simple two-bubble configuration In the context of a fluid description and 
with the emphasis on identifying the basic physics. A simple plecev»is« 
constant density profile is used and the electric field inside and outside 
the bubbles is obtained analytically by the method of images. The method 
is then applied to multi-bubble configurations. In a sense, the present 
work is a non-trivlal extension of the single bubble work of Ossakow and 
Chaturvedl (1978). 

The scope of the present paper is limited to discussing the morphology 
of multi-bubble systems and the time-dependent evolution is not explicitly 
considered. However, the implications of the results will be discussed in 
the context of the behavior of bubbles consistent with the approximations 
used in the analysis. Although we primarily discuss the ESF plasma density 
depletions (bubbles), the technique developed here can be extended 
straightforwardly to the treatment of plasma density enhancements (clouds). 
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In Section II, we develop a two-’-bubble model based on an electrostatic 
(dielectric) analogy and solve the current conservation equation with 
quasi-neutrality for the electric field. Section III describes a multi- 
bubble model and section IV contains the discussion. 

II. A TWO-BUBBLE MODEL 

A. Formulation 

In the present paper, we consider the electric field configuration of 
electrically interacting multi-bubble systems imbedded In a uniform 
background plasma and neutral gas. In order to illustrate the basic 
physics and the theoretical method, we first develop a simple two-bubble 
model. For his purpose, we adopt a sharp-boundary density profile in which 
the plasma density is piecewise constant, being uniform inside (n^) and 
outside (n£) the bubbles and having a discontinuity at the bubble 
boundaries [ Haerendel, 1973; Ossakow and Chaturvedt , 1978; Overman et al. , 
1982]. The bubbles are modelled by two-dimensional cylinders at the same 
altitude with circular cross-sections and the axis of the cylinders arc 
aligned with the earth's magnetic field which is assumed to be uniform 
along the positive z-axis. No neutral wind is included. Figure 1 shows 
schematically the geometry and the coordinate system. The bubbles are 
located at x * - * 0 and x ■ x Q so that the inter-bubble separation distance 
is 2x q . For reference purposes, we denote the cylinders by Cl and C2, 
respectively. The radius of each bubble is a. In the equatorial and low 
latitude regions, the x-axis is along the east-west direction and the 
gravitational force is along the negative y-axis. 

The basic equations in our model are the particle conservation, 
momentum conservation and current conservation equations [see, for example, 
Ossakow, 1981]. In addition, quasi-neutrality can generally be assumed on 
the time scales of interest. In the present paper, we study the morphology 
of two- and multi-bubble systems by solving the current conservation 
equation to obtain the instantaneous electric field perpendicular to the 
magnetic field. 
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The current conservation equation with quasi-neutrality is 

7 • J - 0 (l) 

where J. Is the plasma current density due to Ion and electcon drifts. 
Neglecting the inertia terms in the momentum equations for cold Ion and 
electron fluids, the current density In the lab frame can be expressed as 
(see, for example, Osaakow et al., 1979) 

j + + ( 2> 

where ■ eB/m^c Is the Ion cyclotron frequency, _z Is the unit vector 

along the earth's magnetic field, £ Is the gravitational acceleration, n Is 
the plasma density, and is the lon-neutral collision frequency. In 
(2), we have used the fact that m e /m^ << 1 and have neglected the 

electron j£ x z_ contribution. In addition. In arriving at (2), we have made 
the approximation « !• In the F region, Is typically of the 

order of 10“^ or less. 

The second term In the square brackets In (2) gives the force- 

field (£ and E) aligned drift currents due to the finite lon-neutral 

collisions. It is convenient to separate the electric field £ according to 

~ m i 

£ 3 £ - — t- 

The term (-o^/e) £ Is the component of E cancelling the drift along the 
gravitational field so that the net drift perpendicular to the magnetic 
field Is described by E_. Then, the current J can be written as 


i- °5j 


° 2 V in 


£j = £ + evT- & x £• 
in 









Here, 13 and £ * - g£ are assumed to be uniform. Physically, E^ can be 
thought of as the electric field driving the current J_ perpendicular to the 
magnetic field in the frame moving with the velocity * -8^ £ * z_ 
relative to the lab frame and 0 may be identified as the Pedersen 
conductivity due to v * 0. 

Equation (1) can now be written in the equivalent form 


7 • (nE) - - ^— (& x z) • 7n. (5) 

in 

Perkins et al. (1973) obtained this expression and noted that (5) describes 

a dielectric immersed in a uniform electric field E where 

—o 

So (6) 

in 

In an earlier work, Longmire (1970) utilized a similar magnetostatic 
anology to treat the motion of Isolated ion clouds. Ossakow and Chaturvedl 
(1978) used (5) to analytically study a single bubble system. In this 
dielectric analogy, _E is the polarization (self) electric field of the 
bubbles in the uniform field and E^ corresponds to the total electric 
field (I^j + E) satisfying the boundary conditions across the bubble 
boundaries. 


(oEj), ■ continuous 

(7) 

(Ej)j ” continuous 


and at infinity (x,y ♦ «•) 




♦ e 


( 8 ) 


Note that we have implicitly chosen a reference frame in which the electric 
field of the distant undisturbed ionosphere is Z IJ . The symbols 
I and i refer to the directions parallel and perpendicular to the boundary 
surfaces, respectively. In the present paper, we also adopt the dielectric 
analogy and solve the current conservtlon equation (1) subject to the above 
boundary conditions (7) and (8). For this purpose, it is illuminating to 
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rewrite equation (1) as 


V*(oE.) - 0 
— —u 


(9) 


In the following section, we describe the method of Image dipoles used to 
solve this "dielectric equation". As a matter of notation, Ln the 
remainder of the paper, we use E_ without the subscript J to denote the 
solution of (9). 


B. The Method of Image Dipoles 


The problem of solving Poisson's equation (9) with multiple 
disconnected boundaries Is generally difficult. However, In the case 
treated here with circular cross-sections, the dielectric analogy allows us 
to construct an exact solution. Consider first a single dlelectrLc 
cylinder of radius a centered at x * 0 and a line charge density q located 
at x - b (|b|>a). It is well known ( Smythe, 1968) that the Induced 
electric field outside the cylinder is that of a line charge q'« 

- q(l-K)/(l+K) located at x * a /b and a line charge -q' located at x»0. 
The quantity K is the ratio of the dielectric constant (o^) Inside the 
cylinder to that outsLde the cylinder 


K = 



( 10 ) 


The Induced electric field Inside the cylinder Is that of a single line 
charge q" ■ 2q/(l-HC) located at x - b. If we replace the line charge q by 
a line dipole moment P^ * P^x which Is equivalent to two equal and opposite 
line charges separated by a vanishingly small distance, then we find that 
the electric field due to polarization of the cylinder Is that of a single 
image dipole IP given by 


P , . llZ*) 
- '1+K ' 




(ID 


located at x - a 2 /b. Note that no Image dipole Is present on the axis (to 

be contrasted with the line charge case) and that P and P are collnear. 

—o — 

pointing In the opposite directions. The Induced electric field Inside the 
cylinder Is that due to a dipole moment £* given by 
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located at x * b. 


(12) 


P* ■ p 

- 1+K —o 


We now consider two identical dielectric cylinders of radius a 
centered at x * - x Q and x ■ x Q (see figure l) immersed in a uniform 
electric field E^. Suppose, for the moment, that the two cylinders are 
non-interacting• Then, for the purpose of calculating the polarization 
electric field outside the cylinders, each cylinder may be replaced by a 
dipole moment located at x - -x 0 and x • x Q , where 


1 (ill) a 2 E 

2 '1+K 3 c 


(13) 


and K is defined by equation (10). The components of the self electric 
field are 


E x " 2 p 0 L f ( x+ v y) + f ( x_x 0 » y ^* 


(14) 


and 


E y " 2p 0 [ h < x4x 0 ’ y) + h < x ~ x o' y) J» 


(15) 


where 


f (*.y) 


. x 2 - y 2 
(x 2 +y 2 ) 2 


and 


(16) 


h(x,y) 


2 x y 
(x 2 +y 2 ) 2 


(17) 


Here the x- and y-components of the electric field due to a dipole P Q 
at x - x Q is 2P 0 f(x-x 0 ,y) and 2P Q h(x-x 0 ,y), respectively. The electric 
field E* inside the bubbles is 


E* 


1+K 


(18) 


In the remainder of the paper, asterisks will be used to denote the 
electric fields Inside the bubbles and dipole moments producing the fields. 


7 










We now allow £h^ two cylinders to Interact with each other. In 
addition to the uniform external field E , each cylinder experiences the 

dipole field of the other. As a result, each cylinder Is further polarized 

giving rise to linage dipoles as given by (11) and (12). Iterating the 

method of linage dipoles described above, it Is straightforward to show that 

the total field outside the bubbles Is given by 

E - - E + 1L 2P [f(x+x ,y) + f(x-x ,y)J, (19) 

x o n*0 n L v n ,,7J 

and 

E y “n-0 2p n [ h <x+* n ,y) + h ( x_x n *y)J. < 20 ) 


where f and h are defined by (16) and (17). Here, for n * 0, 



_ . t 1“K\ a 2 _ 

P n = ' W ^ P n-1 * 
n 

(21) 


b n H X o + Vl* 

(22) 

and 

a 2 

X - -- - - . 

n o x +x 

o n-l 

(23) 

For n ■ 0, we have 

x n » x Q and P Q is given by (10). 


Similarly, the total electric field inside each cylinder. 

say, (C2) 

located at x • x Q , 

Is found to be 


E* - 
X 

- TS E o + llit n!o 2P n £(x+ *n- 1,) - 

(24) 

and 

jo 2p „ h <» + «„'>'> • 

(25) 


where P n and x n are defined above. For the other bubble (Cl), the field Is 
obtained by replacing x 0 with -x Q In the functions f and h. 
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It can be seen from (19)-(23) that the convergence properties of these 
series depend on the parameter s defined by 


= il K 

8 = 1+K 


(2x )' 


(26) 


Each successive line dipole moment Is reduced by a factor of s and a 
geometrical factor of order unity from the preceding one. Since x Q > a, we 
have Is| < 1/4 for any K. As a result, the series are generally rapidly 
convergent except for very small center-to-center separation distances 


(x Q ~ a). 


C. Applications 


In order to apply the above dielectric results to the Ionospheric 
bubble and cloud problems In accordance with the dielectric analogy 
(equations (5) and (6)), we identify the dielectric constant o with the 
Pedersen conductivity defined In Section II and recall that E^ is given by 
equation (6). Note that the Pedersen conductivity is proportional to the 
plasma density so that 



where n^ and ^ are the plasma densities inside and outside the bubbles. 
With these Identifications, equations (19), (20), (24) and (25) describe 

the electric field of a two-bubble system In the frame moving with V. - 
-1 

- Si ^ £ x £ relative to the earth. Note that our formalism guarantees that 
the field components obtained above exactly satisfy equation (9) and the 
boundary conditions (7) and (8), as can be verified easily. Thus, the 
solution Is unique. It is also worth noting that all the Image line dipole 
moments are aligned with the x-axls and no higher multipoles such as 
quadrupole moment arise. If the bubbles are at different heights, the 
Interaction can still be expressed as a series of Image dipole moments that 
are parallel to the x-axls but the images are not induced on the x-axis 
Itself. 

As a result of the polarization electric field, the plasma bubbles 
(density depletions) E x » drift. The drift velocity relative to the 
distant undisturbed Ionosphere where the electric field is Eg is given 
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2 

by V - c (E* - E ) x B/B . In particular, the single-bubble rise 
— — —o — 

velocity is (using equation 18) 



For bubbles, K < 1 and is upward. This result has been obtained by 
Osgakow and Chaturvedi [1978]. Because £ is assumed to be uniform, a 
single isolated bubble maintains its circular cross-section as it rises 
with the constant velocity V^. The presence of a second bubble, however, 
modifies the rise velocity significantly. In particular, it is no longer 
uniform and the cross-sections do not remain unchanged. We have 
numerically carried out the summations indicated in (19), (20), (24), and 
(25). The results are depicted in figure 2. The solid lines represent the 
self electric field lines or equivalently the Pedersen current lines while 
the dashed lines correspond to the Instantaneous E x J1 drift 
velocity (Vg) for a two-bubble system (2x Q » 2.5a and K ■ 0 for 100% 
depletion). This figure shows only one quadrant; the actual system is 
spatially uniform in the z-dlrection and is symmetric about the y-z plane 
and the x-z plane. The solid electric field lines are such that the line 
density is proportional to the field strength. This figure clearly shows 
that the electric field inside the bubble is significantly modified from 
the uniform field of an isolated bubble (see (18)). As a general remark, 
the interaction vanishes as K -*■ 1 and the separation distance increases to 
infinity. 

In this example, the electric field strength at the point A (x * x G -a, 
y ■ 0) is approximately 1/3 of that at B (x * x Q +a, ^ *0). The 
prominently nonuniform electric field inside the bubbles has a number of 
Important implications for the dynamic behavior of the two-bubble system. 
The drift lines (dashed) in figure 2 show that different regions of a 
bubble undergo drift in the east-west (x) direction with respect to the 
undisturbed ionosphere. This horizontal drift can be a significant 
fraction of the vertical rise velocity at some points inside the bubble. 
For example, |E*/(E*-E q )| - 0.55 at r - a_, and 6 - ir/8 where r is the 
radial distance from x“x Q , a_ is just inside the bubble surface and 8 is 
measured from the point A (6 - 0). Here, x Q - 1.25 and K » 0 have been 
used (figure 2). The reason for the strong divergence in the field lines 
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2 

x/a 


3 


4 


A drawing of the electric field lines (solid lines) and the x B_ 
drift lines (dashed lines) showing one quadrant of a two-bubble 
system. Tne system is symmetric about the x-z plane and y-z 
plane. The separation distance Is 2x Q ■ 2.5a and K - 0 (100% 
depletion). The points A and B are inside the bubbles at * ■ x Q -a 
and x * x Q +a, respectively. 





Is that the dipole field due to one bubble opposes the internal field of 
the other bubble* 

Figure 3 shows the two-bubble rise velocity (V y ) at a number of points 
inside the bubble (C2) relative to the uniform rise velocity Vj of a single 
bubble given by (27). The ratio R 2 = V y /V^ is plotted for the interior 
points A and B (figure 2) as a function of x Q /a. We note that the rise 
velocity V y is strongly affected by the neighboring bubble for center-to- 
center separation distances (2x Q ) smaller than 5a to 6a. Moreover, 
comparing V y at A and at B, we see that the E x JJ drift velocity is 
sheared. This, together with the conclusions of the preceding paragraph, 
implies that an initially circular cross-section would not remain circular 
so that the two-bubble system described here does not correspond to a 
steady-state system. This should be contrasted with single-bubble models 
in which steady-state solutions are possible for piecewise constant density 
profiles. Because a neighboring bubble generally introduces nonuniformity 
in the electric field, it seems difficult to construct a steady-state two- 
bubble system unless x Q /a is large. 

In interpreting the results, we note that our results do not include 
the time-dependence. The above field configuration exists if the two- 
bubble system as described is created at some time (say, t ■ 0). Thus, it 
would be appropriate as a consistent initial condition for the purpose of 
studying the subsequent time-evolution. In practice, the drift lines shown 
in figure 2 are expected to closely approximate the actual evolution for 
some period of time after t ■ 0 until the distortion changes the field 
topology significantly. However, the general feature of the field with the 
weakest x-component and hence the slowest rise velocity in the region 
nearest to the neighboring bubble should remain unchanged in time. 


III. A MULTI-BUBBLE MODEL 

In the preceding section, we have discussed in detail a two-bubble 
model. In considering a multi-bubble system that may be applicable to the 
ESF phenomenon, the basic physics and the theoretical treatment remain 
unchanged. However, as the number of bubbles Increases, so does the number 
of image dipoles. For an N-bubble system, the n c ^ order expression must 


13 











Pig. 3 The ratio R 2 of the vertical rise velocity V y of a two-bubble 
system to the single-bubble rise velocity (equation (27)) as a 

function of x Q /a. The separation distance Is 2x Q “ 2.5a and K * 0 
(100% depletion). The curve A corresponds to the point A (x ■ x Q - 
a) and the curve B to the point B(x • Xq+s) In figure 2. 
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Include N(N-l) n image dipoles. Fortunately, the Influence of a bubble 
decreases as the inverse square of the separation distance. Thus, the 
nearest and the second nearest neighbors are expected to have the dominant 
effects. A consideration of equations (11) and (12) shows that the second 
nearest neighbors have effects of the order of 4 ^E Q " 0.06 E Q and that the 
third nearest neighbors have effects of the order of 6 E Q * 0.03 E Q for a 
given value of a/x Q < 1. This means that the inter-bubble (and also inter- 
cloud) interaction is short-ranged and that the third nearest neighbors and 
beyond have no significant influence. Thus, only a small number of bubbles 
are necessary to model an N-bubble system with N » 1. Note that the 
influence of the bubbles still vanishes at infinity, differing from systems 
satisfying periodic boundary conditions. 

The theoretical treatment descr.bed is exactly applicable to any N. 
However, in the remainder of this section, we include up to a total of 5 
mutually interacting bubbles. The mathematical manipulations involved are 
analogous to those of the preceding section, resulting in series 
expressions similar to (19), (20), (24), and (25). As before, only dipole 
moments, not higher multipoles, are Induced. Since no new insight is to be 
gained by examining the actual expressions, we give below only the 
results. In figure 4, we shcn the electric field and drift configurations 
of a system with three plasma depletions. The three depletions are again 
modelled by cylinders of radius a, located at x ■ - 2x Q , 0, and 2x Q . Only 
one quadrant is shown. The neighboring bubbles are separated by a distance 
2x Q ■ 2.5a, and are 1002 depleted (K “ 0) as before. The ’’external” 
electric field is E^ given by (6), and the density profile is piecewise 
constant. The solid lines represent the polarization electric 

A 2 

field without E . The quantity c (E - E ) x B /B is then the 
— o — — o — 

instantaneous drift velocity relative to the distant undisturbed ionosphere 
and is represented by the dashed lines. Although the three-bubble system 
is different from a two-bubble system in that the former has a central 
bubble about which the system is symmetric, the general features of the 
field and drift configurations are similar as can be seen by comparing 
figures 2 and 4. That is, the field lines and drift lines inside and 
around the end bubbles of the three-bubble system (figure 4) are similar to 
those of the two-bubble system (figure 2) because the dominant influence 
arises from the nearest neighbors. The central bubble In the three-bubble 









I 


2 

x/a 


3 


A drawing of the electric field lines (solid) and the E. x B. drift 
lines (dashed), showing one quadrant of a three-bubble system* 
The bubbles are placed at x ■ + 2x Q and x ■ 0. The points A, B, 
and C are Inside the bubbles at x - 0, x • a, and x ■ 2x Q +a. The 
separation distance Is 2x Q - 2.5a and K - 0 (100% depletion). 
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system Is affected by two neighboring babbles with comparable effects• The 
field is markedly reduced in the region nearest a neighboring bubble (point 
B in figure 4). 

Figure 4 also shows the distortion in the electric field which renders 
the system non-steady-state, as in the two-bubble case. The drift velocity 
(dashed lines) has an east-west (horizontal) component that may be a 
significant fraction of the vertical velocity. For example, 

|E*/(E*-E q )! • 0.48 at 9 * tt/ 8 from the point B just inside the boundary of 
the central bubble and |E*/(E*—E q ) I * 0.55 at 9 m it/8 from the point x » 
x Q -a and y * 0 just inside the boundary of the side bubble. Thus, in the 
neighborhood of these points, the bubble elements should have significant 
horizontal drifts. 

In addition to the distortion of bubble contours resulting from the 
non-uniform electric field, the figure also shows that the electrLc field 
is substantially reduced from that of an isolated single bubble. This fact 
is Illustrated in figure 5 which gives = I/V^i for the points A, B and 
C corresponding to x =* 0, x * a and x = 2x Q + a, all just InsLde the bubble 
surfaces. Here, V y is the vertical drift velocity of the three-bubble 
system relative to the distant ionosphere. Comparing figure 5 with figure 
3, we note that the qualitative behaviour of the rise velocity is similar 
in both systems, exhibiting significant reduction from that of a single 
bubble system. However, in the three bubble configuration, the electric 
field inside the central bubble is substantially weaker than the two 
neighboring bubbles. Thus, the central bubble has the slowest rise 
velocity. The line D in figure 5 gives the relative vertical rise velocity 
of a five-bubble system calculated at x ■ a_. Curve D, to be compared with 
curve B, shows that the Influence of the additional bubble on the field 
inside the central bubble is small. The reduction in rise velocities is 
increased as the number of bubbles is Increased. However, the influence of 
the bubbles beyond the third nearest neighbor is small. In figure 6, we 
show the relative vertical velocity R * V^/V^ 38 a function of N, the 
number of bubbles, for several values of x Q . For a given separation 
distance, the rise velocity decreases with increasing N and levels off for 
N > 3. Thus, a three-bubble system describes well the basic morphology of 
an N-bubble system. 
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2 3 4 5 

V a 

,g. 5 The ratio Rj of the vertical velocity V y of a three-bubble system 
to the single-bubble rise velocity V^, plotted versus x Q /a. The 
curves A, B, and C correspond to the points A, B, and C In figure 
A. The curve 0 (dashed) corresponds to Rj ■ Vy/V^ at x » a In a 
five-bubble system, to be compared with the curve B. 
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Figures 3 and 3 show that the vertical drift velocity of a bubble Is a 
sensitive function of the separation distance (2x 0 ) except for relatlvly 
large values (2x Q /a >6). In the context of the Raylelgh-Taylor 
instability, this implies that the electric field configuration and the 
drift velocities may depend sensitively on wavelengths (2x Q ). 

As a general remark, we point out that the preceding results derived 
for plasma density depletions for which K < 1 are also applicable to plasma 
density enhancements (clouds) for which K > 1. Calculations for the cloud 
case show that the electric field configurations are qualitatively similar 
to that of the bubble case (figures 2 and 4). In particular, the electric 
field Inside the clouds experiences the greatest reduction in the regions 
facing the neighboring clouds. However, the boundary condition 

K(E ) * (E ) implies that (E. ). is smaller than (E ). by a 

in i out 1 r in i v out L 

factor of for clouds. Thus, for a given separation distance 2x Q , the 
relative distortion in the electric field lines Inside a cloud is less 
pronounced than in a bubble. 

Finally, in figure 7, we have plotted the vertical drift velocity 
versus K 5 n^/^ for a two-bubble (K < l) and two-cloud (K > 1) system. 
The velocity is calculated at x * + (x Q - a) (point A in figure 2) and Is 
upward for bubbles and downward for clouds. The velocity is normalized to 
(equation (27)) for each value of K. Note that the point with K * 1 
does not exist for each line. For K * 1, the ionosphere Is not disturbed 
and there is no bubble or cloud drifting vertically. This is born out by 
the fact that the drift velocity vanishes for a single-bubble (cloud) and 
any multi-bubble (cloud) system. Clearly, the velocities vanish 
differently for different separation distances. Mathematically, the ratio 
of the vertical drift velocity V y - c(E x * - E 0 )/B Q to the single-bubble 
(cloud) drift velocity has the limit 

v^“ “ l 1 “ T?k a 2 f(* + x o ,y>] 

as K + 1. The circled points in figure 7 correspond to the absolute value 
of the quantity which has no physical meaning. 
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{*2 ■ Vy/V^ versus. K for two-bubble (K < 1) and two-cloud (K > l) 
systems. The separation distances (2x Q ) are (a) 2x Q - 2.5a, (b) 
2x 0 ■ 4a, and (c) 2x 0 ■ 10a. 








IV. SUMMARY AND DISCUSSION 

In Che preceding sections, we have solved the current conservation 

equation (1) with quasl-neutrallty for two-bubble and three-bubble systems 
using a dielectric analogy. These two configurations Include the dominant 
near-neighbor Interaction and can model the essential morphology of the 
multi-bubble systems described. This Is demonstrated by actually 
calculating the field Including up to five bubbles. The solutions are 
exact, satisfying the specified boundary conditions on all the multiple 

disconnected boundary surfaces and at infinity. 

Equations (19), (20), (24) and (25) give the solution of the two- 

bubble system as a superposition of Image line dipole moments. Similar 
expressions are obtained for three- and N-bubble (cloud) systems (N > 3). 
In all cases, the interaction is dominated by the nearest neighbors and is 
sensitive to the separation distance (2x Q ) between bubbles. An important 
result Is that the electric field inside the bubbles is generally highly 
non-uniform so that the multi-bubble and multi-cloud systems are not 
steady-state configurations even with piecewise constant density profiles 

(figures 2 and 4). Moreover, the electric field inside the bubbles Is 
significantly weaker than that In a one-bubble system so that the vertical 
drift velocity relative to the undisturbed ionosphere is slower than the 
one-bubble case (equation (27) and figures 3 and 5). For moderately small 
center-to-center separation distances, 2x Q < 6a, the reduction In the rise 
velocity is substantial. This Implies that the electric field 

configuration and the x ji drift velocity may depend sensitively on the 
wavelengths of instabilities causing the Initial density fluctuations 
(e.g., the Raylelgh-Taylor instability). In particular, the bubbles with 
smaller wavelength-to-radius ratio would rise more slowly and be distorted 
more strongly. This conclusion may be particularly applicable to the 
Initial linear or early-time non-linear stages of the evolution. It has 
also been shown that In some regions Inside bubbles, the horizontal drift 
velocity may be comparable to the vertical drift velocity with 

|E*/(E*-E q )| as large as 1/2 (Sections II and III). 

A corollary that follows from the non-uniform electric field in the 
bubbles Is that the polarization Induced J2 x £ drift velocity Is sheared 
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and that a component of the electric field parallel to the local density 
gradient Is developed In the leading edges of the babbles. b'or the simple 
piecewise constant density profiles used for oar analysis, the density 
gradient Is not well defined. In a more realistic density profile, 
however, some previous model calculations indicate that the electric field 
component parallel to the density gradient has stabilizing Influences on 
the Raylelgh-Taylor instability [ Guzdar et al. , 1982] and the £ x li drift 
instability [ Perkins and Poles , 1975; Huba et al. , 1982). lore 
specifically, these calculations show that velocity shear preferentially 
stabilizes the short wavelength modes. In light of these results, we 
suggest that the bifurcation behavior of bubbles and clouds may be 
inhibited by the presence of nearby bubbles (clouds). 

The emphasis of our analysis has primarily been on plasma density 
depletions (bubbles). However, plasma density enhancements (clouds and 
striation fingers) can also be treated in a similar fashion [ Scannapieco 
and Ossakow , 1976; Scannapieco et al. , 1976; Ossakow and Chaturvedl , 
1978]. If a cloud exists in the equatorial or low latitude regie'' , the 
preceding results are all applicable with the replacement of K > 1 (< < 1 
for bubbles) where K - n^/^* In particular, as thi backside an initial 
cloud begins to bifurcate, the small x Q /a and noa 1 I N results may be 
applicable. One point to note is that the electric field inside a two- and 
three-cloud system with a piecewise constant density profile is similar to 
that shown in figure 2 and 4 with weaker field strength in the regions 
facing the neighboring clouds. However, the boundary condition K(E^ n )^ - 
(E 0 u t >i im PH e8 that ( E i n )l I® smaller than (E out ) L by a factor of K~ l for 
clouds. As a result, the distortion in the electric field lines inside a 
cloud is less pronounced than in a bubble. 
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OICY ATTN DOC CON FOR L-389 R. OTT 
OICY ATTN DOC CON FOR L-31 R. HAGER 
OICY ATTN DOC CON FOR L-46 F. SEWAR 

LOS ALAMOS NATIONAL LABORATORY 

P.O. BOX 1663 

LOS ALAMOS, NM 87545 

OLCY ATTN DOC CON FOR I. WOLCOTT 
OICY ATTN DOC CON FOR R.F. TASCHEK 
OICY ATTN DOC CON FOR E. JONES 

OICY ATTN DOC CON FOR J. MALIK 

OICY ATTN DOC CON FOR R. JEFFRIES 

OICY ATTN DOC CON FOR J. 2INN 

OICY ATTN DOC CON FOR P. KEATON 

OICY ATTN DOC CON FOR D. WESTEP.VELT 

SANDIA LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 

OICY ATTN DOC CON FOR W. BROUN 

OICY ATTN DOC CON FOR A. THORNBROUG 

"OICY ATTN DOC CON FOR T. WRIGHT 

OICY ATTN DOC CON FOR D. DAHLGREN 

OICY ATTN DOC CON FOR 3141 
OICY ATTN DOC CON FOR SPACE PROJECT 

SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE, CA 94550 

OICY ATTN DOC CON FOR B. MURPHEY 

OICY ATTN DOC CON FOR T. COOK 

OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON DR. YO SONG 

OTHER GOVERNMENT 

DEPARTMENT OF COMMERCE 
NATIONAL BUREAU OP STANDARDS 
WASHINGTON, D.C. 20234 

(ALL CORRESs ATTN SEC OFFICER FOR) 
OICY ATTN R. MOORE 
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INSTITUTE FOX TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & INFO ADMIN 
BOULDER, CO 80303 

OICY ATTN A, JEAN (UNCLASS ONLY) 

OICY ATTN W. UTLAUT 
OICY ATTN D. CROMBIE 
OICY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
OICY ATTN R. GRUBB 
OICY ATTN AERONOMY LAB C. REID 

DEPARTMENT OF DEFENSE CONTRACTORS 

AEROSPACE CORPORATION 
P.0, BOX 92957 


LOS ANGELES, CA 90009 


OICY 

ATTN I. 

GARFUNKEL 

OICY 

ATTN T. 

SALMI 

OICY 

ATTN V. 

JOSEPHSON 

OICY 

ATTN S. 

BOWER 

OICY 

ATTN N. 

STOCKWELL 

OICY 

ATTN D. 

OLSEN 


ANALYTICAL SYSTEMS ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 

OICY ATTN RADIO SCIENCES 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 
OICY ATTN J. WORKMAN 
OICY ATTN C. PRETTIE 

BOEING COMPANY, THE 
P.O. BOX 3707 
SEATTLE, WA 98124 

OICY ATTN G. KEISTER 
OICY ATTN D. MURRAY 
OICY ATTN G. HALL 
OICY ATTN J. KENNEY 

BROWN ENGINEERING COMPANY, INC. 
CUMMINGS RESEARCH PARK 
HUNTSVILLE, AL 35807 

OICY ATTN ROMEO A. DELIBERIS 

CALIFORNIA AT SAN DIEGO, UNIV OF 

P.O. BOX 6049 

SAN DIEGO, CA 92106 


CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
OICY ATTN D.B. COX 
OICY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 
OICY ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
OICY ATTN H. LOGSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M.H. BORTNER SPACE SCI LAB 

CENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSF, NY 13201 
OICY ATTN F. REIBERT 

GENERAL ELECTRIC TECH SERVICES CO., INC 
HMES 

COURT STREET 
SYRYCUSE, NY’ 13201 
OICY ATTN G. MILLMAN 

GENERAL RESEARCH CORPORATION 
SANTA BARBARA DIVISION 
P.O. BOX 6770 
SANTA BARBARA, CA 93111 
OICY ATTN JOHN ISE, JR. 

OICY ATTN JOEL GARBARINO 
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GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER) 
01CY ATTN T.N. DAVIS (UNCLASS ONLY) 
OICY A’TN TECHNICAL LIBRARY 
01CY ATTN NEAL BROWN (UNCLASS ONLY) 

GTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

OICY ATTN MARSHALL CROSS 

HSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

OICY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

OICY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
400 ARMY-NAVY DRIVE 
ARLINGTON, VA 22202 
OICY ATTN J.M. AEIN 
OICY ATTN ERNEST BAUER 
OICY ATTN HANS WOLFARD 
OICY ATTN JOEL BENGSTON 

1NTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 

OICY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 
P.O. BOX 85154 
SAN DIEGO, CA 92138 

OICY ATTN J.L. SPERLING 

JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAURAL, MO 20810 

OICY ATTN DOCUMENT LIBRARIAN 
OICY ATTN THOMAS POTEMRA 
OICY ATTN JOHN DASSOULAS 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN T. MEAGHER 


KAMAN TEMPO-CENTER FOR ADVANCED STUD: 
816 STATE STREET (P.O DRAWER QQ) 

SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN TIM STEPHANS 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 

LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 

LOCKHEED MISSILES & SPICE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES 6 SPACE CO., INC. 
3251 HANOVER STREET 
PALO ALTO, CA 94304 

OICY ATTN MARTIN WALT DEPT 52-12 
OICY ATTN W.L. IMHOF DEPT 52-12 
OICY ATTN RICHARD G. JOHNSON DEPT 
OICY ATTN J.B. CLADIS DEPT 52-12 

LOCKHEED MISSILE (. SPACE CO., INC. 
HUNTSVILLE RESEARCH 6 ENGR. CTR- 
4800 BRADFORD DRIVE 
HUNTSVILLE, AL 35807 
ATTN DALF. H. DIVIS 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, FL 32805 

OICY ATTN R. HEFFNER 

M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN P. WALDRON 
OICY ATTN L. LOUGHLIN 
OICY ATTN D. CLARK 

MCDONNEL DOUCLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
OICY ATTN N. HARRIS 
OICY ATTN J. MOULE 
OICY ATTN GEORGE MROZ 
OICY ATTN W. OLSON 
OICY ATTN R.W. HALPRIN 
OICY ATTN TECHNICAL LIBRARY SKRVIC? 
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MISSION RESEARCH CORPORATION 

735 STATE STREET 

SANTA BARBARA, CA 93101 


OICY 

ATTN 

P. FISCHER 

OICY 

ATTN 

W.F. CREVIER 

OICY 

ATTN 

STEVEN L. GUTSCHE 

OICY 

ATTN 

D. SAPPENFIELD 

OICY 

ATTN 

R. BOGUSCH 

OICY 

ATTN 

R. HENDRICK 

OICY 

ATTN 

RALPH KILB 

OICY 

ATTN 

DAVE SOWLE 

OICY 

ATTN 

F. FAJEN 

OICY 

ATTN 

M. SCHEIBE 

OICY 

ATTN 

CONRAD L. LONGMIRE 

OICY 

ATTN 

WARREN A. SCHLUETER 


MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01CY ATTN JOHN MORGANSTSRN 
01CY ATTN G. HARDING 
01CY ATTN C.E. CALLAHAN 


PHYSICAL DYNAMICS, INC. 

P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

R & D ASSOCIATES 
P.O. BOX 9695 
MARINA DEL REY, CA 90291 
OICY ATTN FORREST GILMORE 
01CY ATTN BRYAN GABBARD 
OICY ATTN WILLIAM B. WRIGHT, JR 
OICY ATTN ROBERT F. LELEVIER 
OICY ATTN WILLIAM J. KARZAS 
OICY ATTN H. ORY 
OICY ATTN C. MACDONALD 
OICY ATTN R. TURCO 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
OICY ATTN CULLEN CRAIN 
OICY ATTN ED BEDROZIAN 


MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
OICY ATTN W. HALL 
OICY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP. 

1456 CLOVERFIELD BLVD. 

SANTA MONICA, CA 90404 
OICY ATTN E.C. FIELD, JR. 

PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

OICY ATTN IONOSPHERIC RESEARCH LAB 

PHOTOMETRICS, INC. 

442 MARRETT ROAD 
LEXINGTON, MA 02173 

OICY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

OICY ATTN E.J. FREMOUW 


RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

OICY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
80 WEST END AVENUE 
NEW YORK, NY 10023 

OICY ATTN VINCE TRAPANI 

SCIENCE APPLICATIONS, INC. 
P.O. BOX 2351 
LA JOLLA, CA 92033 

OICY ATTN LEWIS M. LINSON 
OICY ATTN DANIEL A. HAML1U 
OICY ATTN E. FRIEMAN 
OICY ATTN E.A. STRAKt*. 

OICY ATTN CURTIS A. SMITH 
OICY ATTN JACK MCDOUGALL 

SCIENCE APPLICATIONS, INC 
1710 COODRIDGE DR. 

MCLEAN, VA 22102 
ATTN: J. COCKAYNE 


i 
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SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 


OICY 

ATTN 

DONALD NEILSON 

OICY 

ATTN 

ALAN BURNS 

OICY 

ATTN 

G. 

SMITH 

OICY 

ATTN 

L. 

,L. COBB 

OICY 

ATTN 

DAVID A. JOHNSON 

OICY 

ATTN 

WALTER G. CHESNUT 

OICY 

ATTN 

CHARLES L. RINO 

OICY 

ATTN 

WALTER JAYE 

OICY 

ATTN 

M. 

BARON 

OICY 

ATTN 

RAY L. LEADABRAND 

OICY 

ATTN 

G. 

CARPENTER 

OICY 

ATTN 

G. 

PRICE 

OICY 

ATTN 

J. 

PETERSON 

OICY 

ATTN 

R. 

HAKE, JR. 

OICY 

ATTN 

V. 

GONZALES 

OICY 

ATTN 

D. 

MCDANIEL 

STEWART 

RADIANCE LABORATORY 


UTAH STATE UNIVERSITY 
1 DE ANGELO DRIVE 
BEDFORD, MA 01730 
01CY ATTN J. ULWICK 

TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 

OICY ATTN W.P. BOQUIST 

TRW DEFENSE A SPACE SYS GROUP 
ONE SPACE PARK 
RF.DONDO BEACH, CA 90278 
OICY ATTN R. K. PLEBUCU 
OICY ATTN S. ALTSCHULER 
OICY ATTN D. DEE 

VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
OICY ATTN W. REIDY 
OICY ATTN J. CARPENTER 
OICY ATTN C. HUMPHREY 
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IONOSPHERIC MODELING DISTRIBUTION LIST 
(UNCLASSIFIED ONLY) 

COPY TO EACH OF THE FOLLOWING PEOPLE (UNLESS OTHERWISE 


PLEASE DISTRIBUTE ONE 
NOTED) 

NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 2037S 
DR. P. MANCE - CODE 4101 
DR. E. SZUSZCZEW1CZ - CODE 4187 
DR. J. GOODMAN - CODE 4180 
DR. P. RODRIGUEZ - CODE 4187 

A.F. GEOPHYSICS LABORATORY 
L.G. HANSCOM FIELD 
BEDFORD, MA 01730 
DR. T. ELKINS 
DR. W. SWIDER 
MRS. R. SAGALYN 
DR. J.M. FORBES 
DR. T.J. KENESHEA 
DR. J. AARONS 
DR. H. CARLSON 
DR. J. JASPERSE 


NASA 
Code 961 

GODDARD SPACE FLIGHT CENTER 
CREENBELT, MD 20771 
DR. R.F. BENSON 
DR. K. MAEDA 
Dr. S. Curtis 
Or. M. Dubin 

NATIONAL TECHNICAL INFORMATION CENTER 
CAMERON STATION 
ALEXANDRIA, VA 22314 
12CY ATTN TC 

COMMANDER 

NAVAL AIR SYSTEMS COMMAND 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20360 
DR. T. CZUBA 


CORNELL UNIVERSITY 
ITHACA, NY 14850 
DR. W.E. SWARTZ 
DR. R. SUDAN 
DR. D. FARLEY 
DR. M. KELLEY 

HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MA 02138 
DR. M.B. McELROY 
DR. R. LINDZEN 

INSTITUTE FOR DEFENSE ANALYSIS 
400 ARMY/NAVY DRIVE 
ARLINGTON, VA 22202 
DR. E. BAUER 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
PLASMA FUSION CENTER 
LIBRARY, NW16-262 
CAMBRIDGE, MA 02139 


COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIECO. CA 92152 

MR. R. ROSE - CODE 5321 

NOAA 

DIRECTOR OF SPACE AND ENVIRONMENTAL 
LABORATORY 
BOULDER, CO 80302 
DR. A. CLENN JEAN 
DR. G.W. ADAMS 
DR. D.N. ANDERSON 
DR. K. DAVIES 
DR. R. F. DONNELLY 

OFFICE OF NAVAL RESEARCH 
800 NORTH QUINCY STREET 
ARLINGTON, VA 22217 
DR. G. JOINER 

PENNSYLVANIA STATE UNIVERSITY 
UNIVERSITY PARK, PA 16802 
DR. J.S. NISBET 
DR. P.R. ROHRBAUCH 
DR. L.A. CARPENTER 
DR. M. LEF. 

DR. R. OIVANY 
DR. P. BENNETT 
DR. F. KLEVANS 








PRINCETON UNIVERSITY 
PLASMA PHYSICS LABORATORY 
PRINCETON, NJ 08540 
DR. F. PERKINS 

SCIENCE APPLICATIONS, INC. 
1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 
DR. D.A. HAMLIN 
DR. L. LINSON 
DR. E. FRIEMAN 

STANFORD UNIVERSITY 
STANFORD, CA 94305 
DR. P.M. BANKS 


UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 
DR. N. ZABUSKY 
DR. M. BIONDI 
DR. E. OVERMAN 

UTAH STATE UNIVERSITY 
4TK AND 8TH STREETS 
LOGAN, UTAH 84322 
DR. R. HARRIS 
DR. K. BAKER 
DR. R. SCHUNK 


U.S. ARMY ABERDEEN RESEARCH 
AND DEVELOPMENT CENTER 
BALLISTIC RESEARCH LABORATORY 
ABERDEEN, MD 

DR. J. HEIMERL 


UNIVERSITY OF CALIFORNIA, 
BERKELEY 

BERKELEY, CA 94720 
DR. M. HUDSON 


UNIVERSITY OF CALIFORNIA 
LOS ALAMOS SCIENTIFIC LABORATORY 
J-10, M3-664 
LOS ALAMOS, NM 87545 
M. PONGRATZ 
D. SIMONS 
G. BARASCH 
L. DUNCAN 
P. BERNHARDT 

UNIVERSITY OF CALIFORNIA, 

LOS ANGELES 
405 H1LLGARD AVENUE 
LOS ANGELES, CA 90024 
DR. F.V. CORONITI 
DR. C. KENNEL 
DR. A.Y. WONG 


UNIVERSITY OF MARYLAND 
COLLEGE PARK, MD 20740 
DR. K. PAPADOPOULOS 
DR. E. OTT 
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